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The wavelength at which a dune pattern emerges from a 
flat sand bed is controlled by the sediment transport satu- 
ration length, which is the length needed for the sand flux 
to adapt to a change of wind strength. The influence of the 
wind shear velocity on this saturation length and on the sub- 
sequent dune initial wavelength has remained controversial. 
In this letter, we present direct measurements of the satu- 
ration length performed in a wind tunnel experiment. In 
complement, initial dune wavelengths are measured under 
different wind conditions - in particular after storms. Us- 
ing the linear stability analysis of dune formation, it is then 
possible to deduce the saturation length from field data. 
Both direct and indirect measurements agree that the sat- 
uration length is almost independent of the wind strength. 
This demonstrates that, in contrast with erosion, grain iner- 
tia is the dominant dynamical mechanism limiting sediment 
transport saturation. 

1. Introduction 

Understanding how wind transports sand is important in 
many geomorphological problems including dunes dynam- 
ics. It is well known that for a given wind velocity u, greater 
than a threshold velocity « t h, the sand flux q saturates to a 
steady value q Ba , t , which results from a subtile interaction be- 
tween grains motion and wind velocity. Although the origin 
of the saturation flux is now well understood, the transient 
regime and the way q adjusts to q sa , t still remains contro- 
versial. If a change of velocity takes place, the flux does 
not adjust instantaneously to its equilibrium value, and a 
phase lag exists between q and velocity it„. Whatever the 
dynamical mechanisms responsible for this lag, it has been 
proposed to encode it into a single length -Lsat, called the sat- 
uration length. By definition, L sat is the length over which 
the flux q relaxes to its equilibrium value q sa .t- Around the 
saturated regime, this relaxation can be modeled by a first 
order equation: 



^ B!lt dx = 9sat _ q ' ^ 

Different physical mechanisms can be responsible for this 
saturation length e.g., the grain hop length [Charm 2006], 
the length needed to accelerate new grains (called the drag 
length) [Andreotti et al. 2002, Hersen et al. 2002, Andreotti 
and Claudin 2007], the length needed to expel new grains 
from the sand bed [Sauermann et al. 2001 ; Parteli et al. 
2007], the length needed for the negative feedback of trans- 
port on the wind to take place [Andreotti, 2004]- Transport 
saturation is limited by the slowest of these processes. 



Copyright 2008 by the American Geophysical Union. 
0094-8276/08/$5.00 



Understanding the physical origin of sand flux saturation 
is crucial to better describe the formation of dunes. Dunes 
result from the interaction between the wind and their shape 
[Bagnold 1941, Lancaster et al. 1996, Wiggs et al. 1996, 
McKenna-Neurnan et al. 1997, Andreotti et al. 2002]. Re- 
cent models have shown that the initial formation of dunes 
is controlled by the balance between three different mech- 
anisms. The first one is the phase-lag between the dune 
elevation profile h and the basal shear velocity u, [Jackson 
and Hunt 1975]. This purely hydrodynamical mechanism, 
related to fluid inertia and dissipation, is the motor o the 
instability and tends to destabilise the bed [Richards 1980, 
Andreotti et al. 2002, Elbelrhiti et al. 2005, Fourriere et 
al. 2008]. On the other hand, two stabilizing effects exist. 
The dynamics of sediment transport tends to stabilise short 
wavelengths, due to the space-lag between sand flux q and 
shear velocity u* . Gravity has a further stabilizing effect, re- 
lated to the dependence of the threshold shear velocity Wth 
on the local slope [Rasmussen et al. 1996]. A stability anal- 
ysis taken into account the three ingredients, [Andreotti and 
Claudin 2007] shows that the most unstable wavelength is 
proportional to L sat and decreases with u*/u th . Saturation 
length then controls the initial size of dunes. 

In this letter, we aim to investigate the physical origin 
of the sand transport saturation length by studying its de- 
pendence on the wind strength. Two methods are used. A 
direct one by measuring the saturation of the sand flux in 
a wind tunnel, and an indirect method based on field mea- 
surements of the instability of sand bed. 
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Figure 1. Spatial variation of the sediment flux over 
a flat sand bed, measured in a wind tunnel for differ- 
ent wind speeds: u* = 0.25m/s (o), u* = 0.33m/s (□), 
u* = 0.34m/s (A) and u* — 0.42m/s (o). In experiments 
(A) and (o) a small input flux was injected to initiate 
saltation. Solid lines: exponential fit (see text). Field 
measurement Elbelrhiti et al. 2005 for d — 185/xm and 
u, ~ 1.3-ttth are also plotted (•). 
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Figure 2. (a) Barchan instability in the Atlantic Sahara due to a storm characterised by a strong dry wind coming from 
inland (Chergui). Aerial photograph taken on january the 5 th , 2005 - credits DigitalGlobe. (b) Wind direction measured 
in Tan- Tan between december the 1 st , 2004, and january the 10 th , 2005 - regular trade winds are along 6 ~ 23°. (b) 
Corresponding values of the basal shear velocity w* . (d) Aerial picture of a transverse dune field in Northwestern Ar-Rub- 
Alkhali (S. Arabia), (e) 220m long dunes in Northwestern Ar-Rub-Alkhali (S. Arabia) (typical local wind: w* ~ 1.2u t h). 
(f) 30m long dunes in the Atlantic Sahara (typical local wind: it* ~ 1.8u t h). 



2. Wind tunnel measurements 

We first report direct measurements of L sa t obtained 
from controlled experiments. They have been performed 
in the wind tunnel of the Cemagref in Grenoble (lm wide, 
0.5m high, and 4.5m long), using sand grains of diameter 
d — 120±40/im from the Hostun quarry [Cierco et al. 2008]. 
An initial 3cm thick sand bed is prepared and flattened with 
a moving bar. This thickness is gently matched to the wind 
tunnel rough bottom over the first ten centimetres. The bed 
elevation profile h(x, t) is measured at regular time interval 
with a vertical resolution of 500/im, turning off the flow. Us- 
ing the conservation of mass, the sediment flux q at position 
x is simply deduced from the erosion data as: 

q(x) = ~ f X d t h{i) di (2) 

JO 



The curves q(x) obtained for different wind strengths are 
depicted in figure 1. In all the cases, the sediment transport 
increases downstream, starting from a null or small input 
flux, and further saturates to a value g sa t, which increases 
with the wind strength. The evolution of q can be divided 
in two phases. A first exponential increase followed by a 
relaxation phase toward the saturation. The initial phase is 
linked to ejection of grains, each saltating grain ejects few 
other grains when it collides the bed. This results into an ex- 
ponential increase of the flux. This process is noticeable for 
weak wind, and when the flux is low [Andreotti and Claudvn 
2007]. This regime is a priori not taken into account in eq. 
1, which aims to describe the relaxation close to the satu- 
ration level. To determine the saturation length, we have 
thus analysed the zone where the flux is larger than half its 
saturated value g sat . The solid lines in figure 1 show the best 
fit by an exponential law of the form g sat [l — e~ ( - x ~ x °^ Lsit ], 
which is the solution of equation (1). As expected, g S at in- 
creases with wind strength. In contrast, the fitted values 
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Figure 3. Measured most unstable wavelength A as a 
function of the rescaled wind shear velocity u*/uth, for 
a mean grain size d ~ 185/im. The star symbols corre- 
spond to superimposed structures on barchans generated 
by stormy events. The circles correspond to low ampli- 
tude transverse dunes generated by averaged winds. The 
solid line is a fit by the model derived in Fourriere et al. 
2008 assuming that the saturation length is constant - 
it gives L sat ~ 68cm. 



of L sat , on the order of 0.7m, are remarkably insensitive to 
those of u* , as well as to the entrance conditions (Fig. 4) . 

Figure 1 also displays a flux profile measured in the field 
by Elbelrhiti et al. 2005 in the Atlantic Sahara. A 20m 
long sand bed composed by grains of size d = 185/im was 
prepared with a bulldozer flatten. The erosion was mea- 
sured after 24 hours of wind fluctuating around the trans- 
port threshold. To determine a characteristic wind strength 
in such a case, we start from velocity time series U(t) mea- 
sured at some altitude. As shown by Ungar and Haff 1987 
and Andreotti 2004, the saturated flux q B!lt is proportional to 
(U 2 — U? h ) H(U — (7th), where Uth is the threshold value for 
transport, and where TL{x) is the so-called Heaviside func- 
tion [H(x > 0) = 1 and H(x < 0) = 0]. We then define the 
effective velocity as that which would give the same mean 
flux, if the wind was not fluctuating: 



«th 



= 1 + 



u(t)\- 
u th J 



H(U(t) - U th ) , (3) 



where the (...) denotes an average over time. Interestingly, 
the ratio so defined does not depend on the height above 
the soil at which the anemometer is placed, which justifies 
that it is also the ratio of the typical shear velocity u* to 
its threshold value wth- For this particular measurement, 
we find if* ~ 1.3u t h- By fitting the curve q(x) in the zone 
where q > q sat /2, we find a saturation length of 1 ± 0.2m, 
twice smaller than that obtained by fitting the whole profile 
(1.7m). The saturation lengths measured here are also sig- 
nificantly smaller than that (L B at ~ 2.3m) which can be de- 
termined from the measurements of Bagnold 1941 (Fig. 62, 
p. 182) with grains of diameter d = 240^im. This is prob- 
ably due to the fact that the erosion rate was measured by 
means of spring balances below sections of the tunnel, of the 
size of the actual saturation length. In contrast, the space 
resolution of 10cm used in our wind tunnel measurements 
allows to reduce the error bars on L sat to 10%. 
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Figure 4. Saturation length L sa t , rescaled by the inertial 
length p 3 / pf d, as a function of the wind shear velocity 
u*, rescaled by the threshold w t h- Direct measurements, 
performed in a wind tunnel (♦) and in the field (A), 
are compared to those determined from the initial dune 
wavelength (storms: * and slipfaceless dunes (o)) 



3. Initial dune wavelength 

The analysis of wavelength at which dunes form from a 
flat sand bed provides an indirect way to measure the satura- 
tion length. This inverse problem is made possible thanks to 
the increased precision of hydrodynamical calculations giv- 
ing the turbulent velocity field around obstacles of small am- 
plitude [Jackson and Hunt 1975, Richards 1980, Fourriere et 
al. 2008]. Given that the instability requires typically 100m 
to develop, wind tunnels cannot be used to form dunes un- 
der controlled conditions. Fortunately, the surface of large 
dunes precisely behave as flat areas of sand and are thus 
submitted to the primary linear instability [Elbelrhiti et al. 
2005] . The main difficulty is to assess the wind velocity to be 
associated to the formation of given superimposed bedforms. 
Figure 2(a) shows a situation in which this can be achieved 
rigourously. Indeed, the periodic superimposed bedforms of 
wavelength A on the stoss slope of this barchan dune is trans- 
verse to the dune itself. They are formed by a strong storm 
almost perpendicular to the regular trade winds (6 = 23°) 
blowing over the Atlantic Sahara (Fig. 2b). This hot and 
dry anomalous wind, called 'Chergui', is due to the motion 
of the Azores anticyclone over Europe and generates dust 
jets over the ocean. Figures 2(b) and (c) show the two suc- 
cessive storms of this type, both characterised by a peak of 
wind velocity and a rotation of the wind direction. We have 
analyzed several of such events, which occur between 5 and 
10 times per year and are presumably responsible for a part 
of the dynamics of barchan fields [Hersen et al. 2004, An- 
dreotti and Claudin 2007]. Importantly, we were present in 
the field during some of these storms and have been able to 
observe the resulting bed instability. In order to determine 
the most unstable wavelength, we have averaged the spac- 
ing between secondary bed-forms over typically 20 barchan 
dunes. The wind velocity is consistently determined on the 
crest of dunes, taking into account the speed-up factor ~ 1.4 
[Andreotti et al. 2002] with respect to the wind strength 
over the surrounding flat ground (Fig. 2c). It is averaged 
over the period of time during which the wind blows in the 
direction perpendicular to the new born crests, following 
equation (3). This procedure allows to obtain data (stars 
in fig. 3) far above the transport threshold (from 1.5wth 
to 4w t h), since the selected events should be strong enough 
to destabilise the surface of dunes within few days. In this 
range of wind velocity, the most unstable wavelength A max 
is roughly independent of u, . 
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Seeking for situations closer to the threshold of transport, 
we have investigated the case of dunes that present evanes- 
cent avalanche slip faces over long period of time. Figure 2 
shows transverse dunes of this type, in northwestern Ar- 
Rub-Alkhali (S. Arabia). They are composed of grains of 
size d ~ 190pim [Abolkhair 1986], i.e. comparable to those 
in the Atlantic Sahara. We hypothesise that these small 
amplitude dunes remain at the wavelength at which a flat 
sand bed destabilises. They can be considered as submitted 
to the average wind at the scale of the year, as defined by 
equation (3). One finds a mean wavelength A = 130m, in Ar- 
Rub-Alkhali, while the same pattern in the Atlantic Sahara 
rather presents a smaller wavelength A = 35m. Analysing 
time-series of the velocity measured in the airports sur- 
rounding these dune fields, it turns out that the charac- 
teristic wind velocity is significantly smaller in the former 
(it* ~ 1.2u t h) than in the later («» ~ 1.4u t h) location. Be- 
cause of the assumption made, these two points (circles in 
Fig. 3) are less reliable than those measured after the desta- 
bilisation induced by a well identified wind event. Still, two 
independent observations make us confident that the large 
wavelength observed in northwestern Ar-Rub-Alkhali is not 
due to a pattern coarsening process. First, we can not see 
there any superimposed bedform. Second, in the border of 
the transverse dune field, one can evidence isolated dunes in 
an intermediate state between slipfaceless dome dunes and 
well developed barchans. This is clearly the signature of 
a dune close to the minimal length of dunes [Hersen et al. 
2002, Kroy et al. 2002, Parteli et al. 2007, Andreotti and 
Claudin 2007]. In Ar-Rub-Alkhali, they can be as large as 
220m, which is ~ 7 times larger than in the Atlantic Sahara 
(compare panels (e) and (f) of Fig. 2). This confirms the 
sharp increase of the most unstable wavelength close to the 
transport threshold (Fig. 3). Based on the comprehensive 
linear stability analysis derived in Fourriere et al. 2,008, 
which includes the effects of size X/d as well as wind w»/uth 
ratios, we have solved the inverse problem and converted the 
data Amax(^*/wth) into a relation between L sat and u*/u t h- 
Figure 4 shows that this independent determination of L aat 
is in fair agreement with the direct one. 

4. Conclusion 

In this letter we have presented the first accurate mea- 
surements of the sand transport saturation length L sat . 
Both direct and indirect methods consistently show that 
L sat is mostly independent on the wind strength. This result 
points to a very simple interpretation related to grain iner- 
tia. Each grain dislodged from the bed needs some length to 
reach its asymptotic velocity. A crude determination of this 
length can be achieved by solving the equation of motion for 
the horizontal grain velocity component v = 

Ef] f d i t = P - IJ f^-v\(U-v), (4) 

where p s is the density of the grain, pf the density of the 
fluid, U the fluid velocity and Cd the turbulent drag coef- 
ficient. This equation can be rewritten in a dimensionlcss 
form, changing variable from t to the grain position x: 

= |1 - «|(1 - «), with v = £ & x= ^^^x (5) 
dx U ip 3 d 

showing that there is a single length scale in the problem, 
which does not depend on the wind velocity U. As a conse- 
quence, the solution of this equation, x = 5/(1 — 5)+ln(l— v) 
shows a spatial transient over a length proportional to 
p s / Pf d. As this is precisely the scaling law followed by L sat 



[Claudin and Andreotti 2006], the saturation transient of the 
sediment flux is certainly controlled by the grain inertia. 
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